Linearized gravitational waves near space-like 
and null infinity 
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Abstract Linear perturbations on Minkowski space are used to probe numerically 
the remote region of an asymptotically flat space-time close to spatial infinity. The 
study is undertaken within the framework of Friedrich's conformal field equations 
and the corresponding conformal representation of spatial infinity as a cylinder The 
system under consideration is the (linear) zero-rest-mass equation for a spin-2 field. 
The spherical symmetry of the underlying background is used to decompose the field 
into separate non-interacting multipoles. It is demonstrated that it is possible to reach 
null-infinity from initial data on an asymptotically Euclidean hyper-surface and that 
the physically important radiation field can be extracted accurately on J'^ . 



1 Introduction 



Asymptotically simple space-times as defined by Penrose [21] are distinguished by 
their qualitatively completely different structure at infinity. Depending on the sign 
of the cosmological constant X their conformal boundary is time-like, space-like 
or null. These space-times have been studied extensively from the point of view of 
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the initial (boundary) value problem and a lot is known about their properties. In 
particular, space-times which are asymptotically de Sitter (X > 0) are known to exist 
globally for small initial data. The same is true for asymptotically flat space-times 
(A = 0) while for asymptotically anti-de Sitter space-times we only have short-time 
existence results for the initial boundary value problem. On the contrary, initiated 
by an important paper by Bizoh and Rostrowroski [2] and based on numerical and 
perturbative methods, there are now strong hints that these space-times are in fact 
non-linearly unstable (for more details on this topic see the contribution by O. Dias 
in this volume). 

In the present paper we want to focus on asymptotically flat space-times. As we 
already mentioned these space-times exist globally for small enough initial data. 
However, it is still not completely understood in detail, how to characterize the 
admissible initial data from a geometrical or physical point of view. The reason lies 
in the fact that in every space-time with a non-vanishing ADM-mass spatial infinity 
is necessarily singular for the 4-dimensional conformal structure. 

In a seminal paper [10], Friedrich shows how to set up a particular gauge near 
space-like infinity which exhibits explicitly the structure of the space-time near space- 
hke and null-infinity. The so called 'conformal Gaufi gauge' is based entirely on the 
conformal structure of the space-time fixing simultaneously a coordinate system, 
an orthogonal frame and a general Weyl connection compatible with the conformal 
structure. The detailed description of this gauge is beyond the scope of this paper and 
we need to refer the reader to the existing literature [8-1 0, 1 6] . Suffice it to say that the 
conformal GauB gauge is based on a congruence of time-like conformal geodesies 
emanating orthogonally from an irutial space-like hyper- surface. The coordinate 
system, tetrad and Weyl connection are defined initially on that hyper-surface and 
are dragged along the congruence of conformal geodesies so that they are defined 
everywhere. 

It turns out that in this gauge the space-time region near space-hke infinity has a 
boundary consisting of future and past null-infinity together with a 3-dimensional 
'cylindrical' hyper-surface / connecting them. In a certain sense, in Minkowski space- 
time, this cyhnder is a blow-up of the (regular) point f. The 'general conformal 
field equations' (GCFE) express the fact that the conformal class of the space-time 
contains an Einstein metric. They are a set of geometric partial differential equations 
(PDEs) generahzing the standard Einstein equations with cosmological constant. 
When split into evolution and constraint equations within the conformal GauB gauge 
the evolution equations take a particularly simple form and it turns out that the 
cylinder I becomes a total characteristic. This means that the evolution equations 
reduce to an intrinsic system of PDEs on /, i.e., they contain no derivatives transverse 
to /. The intrinsic equations are symmetric hyperbolic on / except for the locations 

: = /rio/^, the 2- spheres where / meets future or past nuU-inflnity. There, the 
equations loose hyperbolicity. It is this feature which is responsible for the singular 
behavior of the conformal structure near space-like infinity. 

In a series of papers [10, 12-14], Friedrich has analyzed the behavior of the fields, 
in particular of the Weyl tensor components and their transverse derivatives, along 
the cylinder /. He has shown that generic initial data lead to singularities at I^. The 
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singularities can be avoided if the initial data satisfy certain conditions, one of them 
being the geometric condition that the Cotton tensor of the induced geometry on 
the initial hyper-surface and all its symmetric trace-free derivatives vanish at the 
intersection of / and the initial hyper- surface. It is still not completely clear what is 
the correct geometric classification of those space-times which satisfy Friedrich's 
conditions (however, see [12, 14, 15]). 

Here, we want to discuss these issues from the numerical point of view. It is 
clear, that the fact that the evolution equations extracted from the GCFE cease to 
be hyperbolic is a troublesome feature in a numerical evolution scheme. We want 
to see how it manifests itself in the simplest scenario we can think of: linearized 
gravitational fields on Minkowski space. If we could not control this situation then 
the use of the GCFE for numerical purposes would not be possible. 

The plan of the paper is as follows: in Sect. 2 and Sect. 3 we give brief summaries 
of the basic analytical and numerical results. In Sect. 4 we discuss ways to overcome 
the singular behavior at /+ while in Sect. 5 we show how we can reach future 
null-infinity. The conventions we use and much of the background can be found 
in [22]. 



2 The spin-2 zero-rest-mass equations 

It is well known [22] that perturbations of the Weyl tensor Cahc^ or, equivalently, 
the Weyl spinor 'Fabcd satisfy the equations for a field with spin 2 and vanishing 
rest-mass: 

y^A'(l>ABCD = 0. (1) 

It is also known [22] that this system of equations suffers from Buchdahl condi- 
tions, algebraic conditions relating the conformal ciu^vature ^abcd of the underlying 
background geometry with the field 

which severely restrict the possible perturbations in any conformally curved space- 
time, rendering (1) inconsistent. 

Therefore, we choose flat Minkowski space-time as our background so that ^abcd 
describes small amplitude gravitational waves propagating in an otherwise empty 
space-time. For a detailed discussion and derivation of the exphcit form of the spin-2 
equations in the present context we refer to [1]. Here, we focus only on the relevant 
points. Starting with the standard Minkowski metric in Cartesian coordinates 

| = T/abdX''dX'' (2) 

where T7ab = diag(l, — 1, — 1, — 1), and performing an inversion at the null-cone of 
the origin 
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X — ; X ' X — ^^bX X 



x-x' ' (X-X) 

puts the metric into the form 



8=^^h^- (3) 



(x-x) 

This metric is singular whenever x-x = 0, i.e., on the null-cone at infinity. So we 
define a conformaUy related metric 



g' = Q^g = T/abdx«dx^ Q = -{x-x) (4) 

which extends smoothly to the null-cone of infinity. Note, that space-Uke infinity 
is represented in this metric as the point x" = 0. In order to exhibit the cylindrical 
structure referred to above we perform a further rescaling of the metric using a 
function Jf(r) = r/i (r), where = (x^)^ -|- (x^)^ -|- (x'')^ and is a smooth fimction 
with jU (0) = 1 . Furthermore, we introduce a new time coordinate by defining x" = 
tK{r). These steps give the final form of the metric 



(5) 



Here, we have denoted the metric on the unit sphere by d(0^. Note, that the metric g 
is spherically symmetric. 

The function determines the 'shape' of in the (f,r) -coordinate system. 
We have chosen the two possibilities /i(r) = 1/(1 +nr) with either n = or n = 1 
with the consequence that J^"^ are represented as either a horizontal (n = 0) or 
diagonal (n = 1) line in a (f , r) diagram, see Fig. 1 . They meet the cylinder / in the 
spheres at (f,r) = (±1,0). We introduce the double null-coordinates (m, v) by 

u=Kt-r, v=ia + r, (6) 

which puts the metric into the form 



1 . . 1 



g = ^dwdv— A^daP'. 



Null-infinity is characterized by the vanishing of one of the null coordinates, the other 
one being non-zero (m = 0, v ^ on and w ^ 0, v = on J'^). Both coordinates 
vanish on the cylinder 7. In Fig. 1 we also display the lines of constant u and v, 
corresponding to radial null geodesies. There are two families of null geodesies, one 
leaving the space-time through J^+, the other entering it through ^~ . 

Since the null coordinates are adapted to the conformal structure, we can see 
clearly that the cylinder is 'invisible' from the point of view of the conformal 
structure. This structure is exactly the same as the representation of space-like 
infinity in the coordinates x^, it appears like a point. However, the differentiable 
structures defined by the (m,v) and the (f,r) coordinates are completely different 
near the boundary r = from the one defined by the x" coordinates. 
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Infinity for fi{r)=l Infinity for ^(r) =1/(1 +r) 




Fig. 1 Tlie neigliborliood of space-like infinity in two different (f,r) -coordinate representations 
corresponding to /x(r) = 1 and ^{r) = I / [l + r) 



In order to derive the spin-2 equations explicitly we introduce a complex null 
tetrad {I" ,n" ,nf ,m") adapted to the spherical symmetry, i.e., 

l"d, = ^ {{l-tK') d, + Kd,) , n"d,, = ^ {{l+tK') d, - Kd,) , 

and m" tangent to the spheres of symmetry. Writing the spin-2 equation in the NP 
formalism, computing the spin-coefficients and finally introducing the 'eth' operator 
S (see [17,20,22]) on the unit sphere puts (1) into the form of eight coupled equations 

(l-fJc'M# + Jc5,.#-(3ff'-(5-fc)M)0^-AiS'#-i, k^l:4, 
{l+tK')d,(j)k^Kd,(j)k + {3K' + {k+l)^)(j)k = liQ(j)k+i, k = 0:3 

for the five complex components of the spin-2 field ^abcDi see e.g. [22]. We exploit 
the spherical symmetry of the background Minkowski space-time even further by 
decomposing the field components (pi;{t,r,9,(j)) into different multipole moments 
using the spin-weighted spherical harmonics jF/m 



l>2-km=-l 



(8) 
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Inserting this expansion into (7) and using the action of S and S' on the spin-weighted 
spherical harmonics [22] the system decouples into a countable family of 1 + 1 
systems indexed by admissible pairs of integers {l,m). We find that each mode </)^'" 
propagates along radial nuU geodesies, the 'inner' modes (^f™, 02™' ^3'" propagate in 
both directions while the 'outer' modes ^q"* ^i"* only propagate along one null 
direction. 

The final step in the preparation of the equations is the spUt into constraint and 
evolution equations, which leads to a system of five evolution equations 

{l+tK')dt(l)o-icdr(l)o = -{3k^ -n)(l)o-H«2(l>i, 

^(<;>2 = ^M«o<^>i-^M«o<fe, (9) 

and three constraint equations 

-2if(9r0i +6r^'0i -2tK'iJ.(j)i + aoiJ.{\ -fK:')02 + 0:2^(1 +fK-')0o = 0, 

-2Kdr(l)2+6rn'(l)2 + aoHil-tK')(l)3 + aoH{l+tK')(l)i=0, (10) 
-2Kdr^3+6rlJ,'^+2tK^IJ,^ + OoHil +tK^)^ + a2M(l -tK^)^4 = 0. 

Note that we have dropped here the superscripts from (^^™ and we introduced the 
quantities Oq = + and a2 = ^/UT+V)^^. Since K-(O) = it is obvious that 

the evolution equations reduce to a system intrinsic to /. Since on / also K"'(0) = 1, it 
follows that the coefficients in front of the time derivatives of 0o (respectively (^4) 
vanish when t = — I (respectively when t = 1). 

In Fig. 2 we show again the neighbourhood of / and ^+ for ju(r) = 1/(1 + r). 
This time we show the characteristics also in the unphysical part of the diagram. 
The non-shaded region bounded partly by the thick broken line is the domain of 
hyperbolicity of the evolution equations, i.e., the domain where f < ±|K''(r)|. Notice 
that every neighborhood of contains regions where the hyperbolicity breaks down. 
Apart from the fact that one carmot hope to get existence and uniqueness of solutions 
beyond that region its presence also makes the numerical evolution challenging. 
For instance, setting up an initial boundary value problem with the left boundary at 
negative values of r does not make sense if the evolution is to reach up to /+. Even 
with the left boundary on / it is not possible to go beyond since the evolution hits 
the non-hyperbolicity region. 
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Characteristics for fi(r) =l/(l +r) 
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Fig. 2 The chai'acteristics of the evolution equations in a neighborhood of /+ in the diagonal 
representation. The shaded region is the domain where the equations fail to be hyperbolic, i.e., 
where t <K'{r). The corresponding neighborhood of is obtained by reflection at the lower border 
of the diagram accompanied by an interchange of u and v. The thick (green) broken line indicates 
the boundary of the domain of non-hyperbolicity. 



3 Numerical methods and tests 

Again, we give a brief summary and refer to [1] for further details. We have different 
equations for different values of the multipole index I. Here, we focus only on the 
case / = 2 without mentioning it any further. We have done evolutions with other 
values of / which lead to very similar results. The evolution equations (9) are solved 
on a spatial interval < r < 1 as an initial boundary value problem. We use the 
method of lines constructing the spatial discretisation using a 4* order accurate finite 
difference scheme. The initial data are obtained either from exact solutions or by 
solving the constraint equations (10) explicitly in terms of two free functions. The 
boundary at r = is a characteristic so we do not need to specify any free functions 
there, while the boundary at r = 1 is an artificial boundary, which needs exactly 
one free function for the component 0o that propagates through this boundary into 
the computational domain. The boundary conditions are implemented using SBP 
operators [6, 18,24] and the SAT penalty method [3,4, 19,23]. The semi-discrete 
system of ODEs is solved using the standard 4* order Runge-Kutta method. 

This code stably propagates the initial data from f = up to any value of f < 1 , 
independently of whether we use the diagonal or the horizontal representation of 
.y . The code converges to 4* order, the constraints propagate and the constraint 
violations remain bounded. However, the two representations behave differently 
when we attempt to reach f = 1 . 

In the diagonal case we can reach t — I exactly and the code converges in 4* order, 
indicating that the numerical problem is still well-posed. With any fixed time-step T 
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we can also step beyond t=ltot=l + T. However, the code fails to converge at 
/ = 1 + T for every T > 0. Clearly, the loss of hyperbolicity at f = 1 is also responsible 
for the lack of well-posedness of the numerical problem. Referring back to Fig. 2 it is 
clear that there will always be a sufficiently high resolution in time and space which 
will detect the domain of non-hyperbolicity of the evolution equations beyond /+. 

In the horizontal case the situation is apparently worse. In this case, it is principally 
impossible to reach t = l simply because the propagation speed of ^4 becomes infinite 
and so stability is lost not only at the point /+ but at the entire time-slice / = 1 which 
coincides with J^+, a characteristic hyper-surface, in this representation. 

However, in principle we can come arbitrarily close to ? = 1 (within numeri- 
cal accuracy) if at every time we choose the next time-step small enough so that 
the CFL criterion — the numerical domain of dependence includes the analytical 
domain of dependence — is satisfied. Such an adaptive time-stepping scheme has 
been implemented and we have in fact demonstrated that we can come very close to 
t = 1 without loosing convergence. However, since the time-steps decrease to zero 
exponentially with the number of steps the simulation takes arbitrarily long. 



4 Beyond /+? 

Given that we can reach ? = 1 in a stable fashion in the diagonal representation we 
may ask the question as to whether it is possible to continue beyond /+? Clearly, for 
the reasons discussed above we cannot simply continue the computation because 
we run into the non-hyperbolicity region near 7+. Referring back to Fig. 2 we see 
that one family of the characteristics — corresponding to the field component ^4 — 
asymptotes to the cylinder and, subsequently, to in a non-uniform way. This 
non-uniformity is due to the fact that the coefficient in front of the time derivative 
of (p4 in the evolution system vanishes at f = 1, causing the loss of hyperbolicity of 
the system. Can we avoid this problem? There are a few possibilities which come to 
mind: 

• Chop the computational domain by dropping points outside from the left. 

• Change the radial coordinate so that becomes the left boundary of the com- 
putational domain. This is a form of ^-freezing (see [7]). 

• Change the time coordinate so that the (space-like) time-slices tilt upwards to- 
wards 0/+. 

We have considered the first two possibilities. These two methods are complementary 
to each other in the sense that in the first case we change the computational domain 
but not the system, while in the second case we change the equations but not the 
computational domain. 

Fig. 2 shows the behavior of the characteristics for <^o and <^4 as well as the 
region in which hyperbolicity fails. The surfaces that we are evolving our data on 
are horizontal and to the right of the vertical solid and broken black lines. Note that 
for f > 1 these surfaces are guaranteed to intersect the region in which hyperbolicity 
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fails. In addition, for surfaces with t > I a new boundary condition for ^4, on the 
'left' of the grid is required. However, since is one of the characteristics for (^4 
this boundary condition cannot influence the physical region to the right of J^+. 

To solve both of these problems we set the values of (|)o, • • • , <|>4 to zero after some 
predetermined grid point that is beyond future null infinity, but before the region 
in which hyperbohcity fails. That is, between the broken green line and J^+. The 
equation for the broken green line is 

and the equation for future nuU infinity, where t> 1, is 

1-^=0. 
r+l 

We chose, after some experimentation, that the 'cut off point, after which all data 
values would be set to zero would be 10% of all grid points beyond That is, the 
values of (^0, ■ ■ ■ ) <?*4 would be set to zero on the 90% of grid points between r = 
and J^^. This condition, of setting values to zero, allows us to numerically provide 
the necessary boundary condition for ^4 and cope with the lack of hyperbohcity in 
the region of the grid to the left of the broken green fine. 

Other than this technicality, the same evolution scheme, with the same fixed 
step size was used to evaluate the initial data up to f = 1.1. We evolved an exact 
solution and estimated the error at ? = 1 . 1 by comparing the values of (^0 ) • • • , <?*4 . on 
the 'physical' portion of the grid to the exact solution. The L2 measure of the error 
is presented in Table 1. We only give the error for the <po and <p4 components. The 
convergence rates for , (j)2 and are, roughly, a linear interpolation between those 
for (j)Q and (p4. Note, that the error in <j)o is of the order of 10~^^ while the error in <j)4 
is roughly 10"^. 



Grid Points 


logzdl^lb) 


Rate 


<^4 

l0g2(l|4||2) 


Rate 


100 


-28.38 




-5.93 




200 


-31.75 


3.37 


-5.98 


0.05 


400 


-35.21 


3.46 


-6.07 


0.10 


800 


-38.71 


3.50 


-8.41 


2.34 



Table 1 Absolute error A compared to an exact solution. norm and convergence rates at time 
t = 1.1 for (j>o,(j>4. The calculation was done with a fixed time-step. 



In the second case we change the r-coordinate in a very simple minded way using 
the coordinate transformation 



rt-^r = r — t + l, 



tt-^t =t, 



fOTt> 1. 
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This has the effect that J^"*" is given as the locus {r = 0}. Note, that the coordinate 
transformation is only continuous and not even Since the partial derivatives 
transform according to 

dt = dt-df, dr = dr 

we still have the problem that the coefficient in front of the time derivative of (^4 
vanishes at f = 1 . Clearly, we cannot make the system regular at f = 1 by a coordinate 
transformation. We can try to avoid the evaluation at f = 1 by arranging the time 
stepping to 'straddle' ? = 1 so that we never actually hit it exactly. However, this 
has the consequence that the 'kink' that we obtain due to the non-smoothness of 
the coordinate transformation induces oscillations at the right boundary. We would 
probably be able to avoid those if we used a smoother coordinate transformation. 
However, this would change the equations in much more complicated ways and we 
have not pursued this any further. 

The third possibility mentioned above — changing the time coordinate — ^has the 
effect that the time-slices globally approach which is exactly the feature that we 
see in the horizontal representation. This prompted us to look at the relationships 
between the different conformal representations in more detail. 

The two representations are related by a conformal rescaling and a coordinate 
transformation. We can derive the corresponding relationships for the field compo- 
nents as follows. Let /i„(r) = 1/(1 +nr) and define & =iii/no. Let go and gi be the 
metrics corresponding to the horizontal and diagonal representations, then we have 

gi = ^go = &ho- (11) 
Mo 

Furthermore, the two time coordinates to and h in the two representations are related 
by 

= rnoto = ruih ^ ti=0 ho- (12) 

The spin-2 field has conformal weight —1 under conformal rescalings. This implies 
that with (11) we also have 

<^1bcd = &-'<^abcd- (13) 

In order to get the behavior of the field components under conformal rescalings we 
observe that (11) impUes that the tetrad vectors rescale as 

1"^=©-HS, etc (14) 

and the spin-frame correspondingly rescales with ©"'Z-^. Taken altogether, these 
transformation properties imply that 

^litur)^©-'<l)^{@-'to,r). (15) 

We can now make use of these relationships in the following way. Suppose we want 
to evolve initial data in the diagonal representation. We rescale the data on the initial 
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hyper-surface ti ~0~to using (15) into initial data for the horizontal representation 
and evolve them with the system for the horizontal representation up to a time 
f 1. Then we undo the rescaling with (15) and obtain the solution in the diagonal 
representation. In Fig. 3 we present the results of these operations. The surface plot in 



(I^ with n=1 with n=1 




too, 0.2 0.4 0.6 0.8 1 



Fig. 3 Comparison between tlie two representations. For detailed explanation see text. 

the upper left shows the component 04 evolved with the diagonal representation n~l, 
while the plot on the lower left shows 04 obtained in the horizontal representation 
using the same data as for the case n = 1 but rescaled according to (15). The contour 
plots on the right show the obtained solutions. Above is the solution directly obtained 
in the diagonal representation while below is the solution obtained after rescaling 
back from the horizontal representation. The contour plots agree visually. Note, that 
the t coordinates in the two surface plots are not the same. They refer to fo in the 
lower plot and to t\ in the other. 

In the horizontal representation we come arbitrarily close to f = 1 . Hence, we can 
'almost' compute the entire space-time (at least in this simple set-up). In this sense 
the horizontal representation is much more efficient than the diagonal one. In Fig. 4 
we show the contour plots of the rescaled solution in the diagonal representation. 
Clearly, it extends way over the f 1 time-slice reaching ^+ to within the graphical 
resolution. 
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Fig. 4 Contour plots of (left) and (j>4 (right) computed with n = but rescaled to the diagonal 
representation. 

5 Reaching null-infinity 

Since the horizontal representation allows us to almost reach null-infinity one may 
wonder whether it might be possible to extend the computation to ^+ from the 
last time-level at f = f , say? To answer this question we refer back to Fig. 1 where 
we showed the space-time in the horizontal representation together with the two 
families of characteristics. The 'outgoing' characteristics (those intersecting 
are well-behaved while the other family asymptotes to / U ^+ in a non-uniform 
way. Referring to Eq. (7) we see that all components except for 04 propagate along 
the well-behaved characteristics. This suggests that we use those four propagation 
equations to estimate the values of (j)Q, 02> 03 on Since we cannot use any 
values of 04 beyond t — twe are forced to use an Euler step to get to f = 1. 

As for 04 we observe that its propagation equation reduces to an intrinsic equation 
on J^+, which we could integrate if we had initial conditions for 04 at r = 0, once 
we know the values of the other components on The equation for 04 on ^+ 
reduces to 

r5,-04 - 204 = a203. (16) 

This equation is singular at r = and the requirement that 04 and its derivative be 
bounded forces the initial condition 

04(1,0) = -y03(l,O). 

So we see that we can perform the last step to J^+, albeit only with a first order 
method. 

We tested this approach by evolving compactly supported initial data obtained 
by explicitly solving the constraints (the same ones as in Sec. 4). The error at f = 1 
was estimated by comparing the values produced in some simulation to the values 
produced in the highest resolution simulation with 6400 grid points. The L2 measure 
of the error is presented in Table 2. Since the 0o and 04 components mark the extreme 
cases, we do not present the error convergence rates for 0i , 02, 03. Their convergence 
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rates are comparable to <j)o. Remarkably, the convergence rate for ^ is close to 4* 





<l>0 




<t>4 




Grid Points 


logzdl^lb) 


Rate 


l0g2(l|4||2) 


Rate 


200 


-7.99 




0.082 




400 


-11.48 


3.49 


-0.44 


0.52 


800 


-14.98 


3.49 


-1.02 


0.57 


1600 


-18.48 


3.50 


-1.73 


0.71 


3200 


-22.07 


3.58 


-2.81 


1.08 



Table 2 Absolute error A compared to a 6400 grid point simulation using the norm and conver- 
gence rates at time f = 1 for ^, (^4 for compactly supported data in the n = space-time. 

order, while for we get the expected first order convergence. The accuracy in the 
radiation component is roughly 10~^ while for (^4 it is not good, only « 10%. This 
low accuracy is due to the singular behavior of the propagation equation (16) for (p4 
at r = 0, which we might not have taken care of appropriately. There is still room for 
improvement because we could use the intrinsic transport equations (see [1,11]) along 
the cylinder — mentioned in the introduction — to drag along as many derivatives of 
(j)4 as we please. These could be used to construct a Taylor approximation for (j)4 near 
r = to get the integration process started more smoothly. 



6 Conclusion 

In this paper we have presented a study of the spin-2 equation in the neighborhood 
of spatial infinity in Minkowski space-time. Since the perturbations of the Weyl 
curvature on flat space obey this equation we can interpret this system as a model 
for small amplitude gravitational waves. We used this model to study the asymptotic 
properties of the fields — and, hence, to some extent also of the perturbed space- 
time — close to spatial infinity. This region is still not completely understood and we 
hope that our work will ultimately contribute to the complete understanding of this 
issue. 

We have shown here that it is possible to generate a complete evolution from 
initial data on an asymptotically Euclidean hyper-surface to the asymptotic regime 
including null-infinity even though the equations show a certain degeneracy at the 
cyhnder which represents spatial infinity in Friedrich's conformal GauB gauge. Using 
the horizontal representation as described in Sects. 4 and 5 we can get the radiation 
field (^0 quite accurately on The drop in convergence for the component (p4 is 
akin to the loss of smoothness of null-infinity in the first result on the global stability 
of Minkowski space by Christodoulou and Klainerman [5]. This was caused by the 
loss of peeling in the corresponding component of the Weyl tensor. 

Of course, we deal here with a highly simplified system and the fact that we 
can make things work here does not immediately imply that it wiU also work in the 
physically relevant 3D cases. However, this toy system does provide valuable insights 
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into the possibilities of and the restrictions imposed by the structure of spatial infinity. 
We should point out that the structure of the fully non-linear general conformal field 
equations is very similar to the linear spin-2 system. This is essentially due to the 
fact that the spin-2 system results from the Bianchi equations obeyed by the rescaled 

Weyl spinor. 

Our next steps will be the removal of the artificial boundary at r = 1. This will 
allow us to perform an entirely global evolution of the mode decomposed spin-2 field. 
The challenge here is to get the centre under control because due to the spherical 
symmetry this will be a singular point for the equations. Then, we intend to remove 
the mode decomposition and look at the linearized system in three spatial dimensions. 
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